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Abstract 

In this review talk I discuss the latest developments in the DGLAP -based global analysis of nuclear parton dis- 
tribution functions (nPDFs), focusing on the recent EPS09, nCTEQ and DSSZ global fits. I also briefly review the 
recent analysis for assigning a spatial dependence to the globally analysed nPDFs, resulting in the new sets EPS09s 
and EKS98s ("s" for spatial). With these, one can now compute nuclear hard-process cross sections and estimate their 
nPDF-originating uncertainties in different centrality classes for the first time consistently with the global nPDF fits. 
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1. Introduction: Global DGLAP analyses of nPDFs 

Inclusive cross sections of hard processes involving a large interaction scale Q » Aqcd in a high-energy hadronic 
or nuclear collision of particles A and B can be computed using the QCD collinear factorization theorem, 

do AB^k+X = y.A ( g 2) 9 y.2J ( g 2) ^ d& U^X' + 0(1/ Q 2 ), (1) 

i,j,X' 

where dcr are the perturbatively calculable partonic pieces, and ff'(f^) is the universal, process-independent, parton 
distribution function (PDF) for a parton flavor i(j) in A(B). The PDFs are of nonperturbative origin, but their evolution 
in the scale Q 2 can be obtained from the DGLAP equations JT] |2j [3] |4] derived from perturbative QCD. 

A precise knowledge of the PDFs is vital for interpreting any hard-process results in p+p, p+A, d+A, and A+A col- 
lisions at the present BNL-RHIC and CERN-LHC colliders. Consequently, global analyses which exploit a multitude 
of experimental hard-process data and the DGLAP evolution, have been developed to determine the nonperturba- 
tive input in the PDFs. Traditionally, in these analyses the PDFs are parametrized at some perturbative initial scale, 
Q - Qo Aqcd, where the power corrections, 0(1 /Q 2 ) in Eq. |T|), can be expected to be small. Performing the 
DGLAP evolution to higher scales, and computing the cross sections of the measured observables (with each param- 
eter candidate-set at the time), one then iteratively determines the parameter set giving the best fit with the measured 
data. A statistical error analysis is performed after this, to determine how the experimental statistical and systematic 
errors translate into uncertainties of the PDFs. With the resulting error sets and best fit for the PDFs, one may then 
compute how the uncertainties in the PDFs propagate into the hard-process cross sections. Excellent fits for the free- 
proton PDFs have been obtained in this way, resulting in the sets like CT10 |5| and MSTW [6|. Also neural network 
techniques have been successfully developed, resulting in e.g. the set NNPDF2.0 Q. 

Challenges for the global PDF analyses are posed, e.g., by the treatment of experimental errors and weighting of 
those data sets which contain a small number of datapoints but which still offer very valuable constraints. Due to the 
large dimension of the fit-parameter space, <9(15 - 30), and due to the complexity of the NLO (and NNLO) cross 
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sections, very fast solvers for the DGLAP evolution and cross sections are also needed. For the global analysis of 
nuclear PDFs, the main further challenge is essentially the fact that we have fewer types of data and less data points 
covering a much more limited x, Q 2 region at our disposal, i.e. less constraints available than in the free proton case. 
In addition, one needs to account for the A dependence, and also the spatial (impact parameter) dependence of the 
nPDFs. All the global analyses so far, listed in Table 1, have been performed for the spatially averaged nPDFs probed 
in minimum-bias nuclear collisions with no cuts on the collision centrality (impact parameter). An obvious weakness 
of these nPDFs then is that it has not been possible to consistently compute nuclear hard-process cross-sections in 
different centrality classes. 



year 


set 


Ref. 


order 


data types used 


error analysis 


1998 


EKS98 


E] ED 


LO 


l+A DIS, p+A DY 


no 


2001 


HKM 


ED 


LO 


l+A DIS 


yes 


2004 


HKN04 


CIl 


LO 


l+A DIS, p+A DY 


yes 


2004 


nDS 


ma 


NLO 


l+A DIS, p+A DY 


no 


2007 


EKPS 


mo 


LO 


l+A DIS, p+A DY 


yes 


2007 


HKN07 




NLO 


l+A DIS, p+A DY 


yes 


2008 


EPS08 


na 


LO 


l+A DIS, p+A DY, 11"=, tt°, tt* in d+Au 


no 


2009 


EPS09 




NLO 


l+A DIS, p+A DY, n° in d+Au 


yes, — > error sets 


2009 


nCTEQ 


lfT71fT8l 


NLO 


l+A DIS, p+A DY 


yes 


2010 


nCTEQ 


m 


NLO 


l+A and v+A DIS, p+A DY 


yes 


2012 


DSSZ 


ED 


NLO 


l+A and v+A DIS, p+A DY, 





7T°, tt* in d+Au, computed with nFFs yes, — > error sets 



Table 1: The developments in the global DGLAP analysis of nPDFs since 1998. The new elements are indicated with underlining. DIS = deep 
inelastic scattering; DY = Drell-Yan dilepton production; nFFS = nuclear fragmentation functions; bold face = reviewed this talk. 

Table 1 summarizes the key developments in the global analysis of the nPDFs. Next, I will briefly review the most 
recent ones of these. The impact-parameter dependent nPDFs of Ref. [21 1 (see also Helenius [[22] in these Proc.) will 
be discussed after this. I will not address the possible power corrections in the cross sections 11231 . or nonlinearities 
in the scale evolution, gluon saturation phenomena or other evolution equations than DGLAP here. It is also good to 
keep in mind that the global DGLAP analyses of the nPDFs do not address the interesting question of the origin of 
the nuclear effects. For such studies, see e.g. the modeling in Refs. [24, 25 , 26 1 and the reviews in Refs. [27, 28 1. 

2. EPS09 

Following the EKS98 framework ]8l |9), in the EPS09 analysis lfl6l the (spatially independent) bound-proton 
nPDFs are defined in terms of the nuclear modifications, Rf(x, Q 2 ), and free-proton PDFs, f?(x, Q 2 ): 

ff\x, Q 2 ) = Rf(x, Q 2 )f t p (x, Q 2 ). (2) 

The bound-neutron PDFs are obtained by assuming isospin symmetry. As the free proton baseline, EPS09 uses the 
CTEQ6.1M free proton NLO PDFs BU). Accordingly, the initial scale in EPS09 is set to Q = 1.3 GeV, and the 
MSbar regularization scheme and zero-mass variable flavour-number heavy-quark scheme (i.e. heavy quarks are 
taken massless and generated only radiatively through the DGLAP evolution above the mass thresholds) are adopted. 
Uncertainties in the free-proton PDFs are not considered. As the constraints, EPS09 exploits 25 data sets from l+A DIS 
[E139 and NMC experiments], 6 sets from DY in p+A [E772 and E866], and one set from inclusive minimum-bias n 
production in d+Au collisions at RHIC [PHENIX], altogether 929 datapoints (for the Refs., see lfl6l ). The quality of 
the EPS09 best fit is excellent: /datapoints = 0.79. The statistical error analysis is done through the Hessian method, 
by first diagonalizing the Hessian matrix and finding the parameter eigendirections as linear combinations of the 15 fit 
parameters. The 15 + 15 error sets, pairwise assigned to each of the uncorrelated eigendirections, are determined with 
a "90% confidence criterion" (see the details in [ 16 1) which in the EPS09 case is equivalent to allowing the x 2 to vary 
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Figure 1 : (Color on-line) Left: The nuclear effects to bound-proton NLO PDFs in a lead nucleus as a function of .v, at two different scales, according 
to the EPS09 global analysis. Comparison with the earlier NLO nPDFs is also shown. The error bands are computed from the 15+15 error sets 
obtained through a "90% confidence criterion". From [16]. Right: Q 2 evolution of the EPS09 nuclear effects at fixed values of x. From [29]. 
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Figure 2: (Color on-line) The EPS09 fit to some of the l+A DIS jF^(x, Q 2 )l\F^(x, Q 2 ) data (upper left), p+A DY 
(jda- pA / clx2<lQ 2 )/ (I dcr^ /dxidQ 2 ) data (upper right) sets used, as well as to the minimum-bias PHENIX nuclear modification factor data for 
7T° production in d+Au collisions, ^(da^^ / dpjdy) l(d(T w I dpjdy), in mid-rapidity at RHIC 1311 (lower panel). The n° production is computed 
with the INCNLO code (32J and NLO KKP vacuum fragmentation functions (FFs) HTI . From fT6l . 

by Ax 2 - 50, i.e. ca. 5%. The resulting error band and comparison with the earlier NLO nPDFs is shown in Fig.[T]for 
the lead nucleus. The good simultaneous fit to the DIS, DY and n {) data is demonstrated in Fig. [2] As seen in Fig.[T[ 
large uncertainties remain especially in the small-x and large-x gluons. The DGLAP evolution diminishes the small-x 
gluon shadowing and also its error band, as seen in the lower panel. The decisive role of the n° data 11311 in fixing 
the mid-x EMC -effect of gluons is demonstrated in the lower panel of Fig. [2] As a summary of the EPS09 analysis, 
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we can conclude that the excellent simultaneous fairly tensionless fit to the different types of nuclear hard-process 
data suggests that factorization works well in the energy range studied, and that the extracted nPDFs seem universal 
in the region x > 0.005, Q > 1.3 GeV. Also the procedure for propagating the experimental errors into the nPDF 
uncertainties and error sets is seen to work well. 

3. nCTEQ 

The latest nCTEQ global analysis for nPDFs, presented in Ref. |fT9ll , is built on a similar set-up as the CTEQ6M 
free proton PDFs. Instead of the nuclear modifications to the PDFs, nCTEQ parametrizes directly the initial bound- 
proton PDFs, jf (x, Qq), at a the scale Qo - 1.3 GeV. Heavy -quark effects are included here in a more sophisticated 
way, using the general-mass variable-flavour number scheme (GM-VFNS). Altogether 31 charged-lepton+A DIS data 
sets [E139, NMC, EMC, E665, BCDMS, HERMES] and 6 DY p+A data sets [E772, E886], consisting of altogether 
708 datapoints were used as constraints to the nPDFs. As a new element relative to EPS09, nCTEQ includes also 8 
v(v)+A DIS data sets [CHORUS, NuTeV, CCFR] with as many as 3134 datapoints into the global analysis. The RHIC 
data have not, however, been used here at all. For all the data Refs., see fl9l . 

The interesting main observation from the nCTEQ analysis is essentially shown by Table 2: with increasing 
weight w to the neutrino-DIS data, the Z ± -DIS&DY fit and global fit deteriorate quite significantly, which indicates a 
tension between the / ± -DIS&DY data and the v-DIS data. We can also see that the fit to the v-DIS data alone is of a 
poorer quality than that to the / ± -DIS&DY data alone. Based on the detailed statistical analysis, the nCTEQ results 
in (19 1 would seem to suggest that the nPDFs are not universal but process-dependent. As pointed out in |[T9l l34l . 
this conclusion depends, however, on the analysis details, in particular on the extent to which the systematic errors of 
the v-DIS data are regarded correlated. Namely, as pointed out in ||35l , where these errors were taken uncorrelated, 
no such a strong conclusion was found. In addition, in Ref. 1 35 1 it was demonstrated that while a relatively good 
agreement with the CHORUS and CDHSW data was found, the NuTeV data, and especially the fluctuations in the 
normalization of the data from one neutrino energy to another, causes the worsening of the fits. This tension problem 
with the NuTeV data alone (which form ca. 2/3 of all the v-DIS data) was also acknowledged in ll34l . At this point, my 
conclusion is that factorization seems to work for the nPDFs, even with the v-DIS data included. To more stringently 
test the universality of the nPDFs, obviously more precise v-DIS data would be needed in the future. 



w 


Z ± -DIS & DY 


v ± -DIS 


total ;^ 2 /datapt 





0.9 




0.90 


1/2 


0.96 


1.40 


1.32 


1 


1.04 


1.36 


1.30 


CO 




1.33 


1.33 



Table 2: The worsening of the global fit as a function of the weight w to the neutrino-DIS data in the nCTEQ analysis. From Ref. 1 19 |. 

4. DSSZ 

The latest global analysis for the nPDFs is presented by DSSZ in Ref. 12011 . Here the baseline free-proton PDF 
set is MSTW 0, and the nuclear modifications Rf{x, Qq) of the bound protons are parametrized at an initial scale 
Qo = 1 GeV. The MSbar regularization and GM-VFNS heavy-quark scheme are used. As constraints for the nPDFs, 
DSSZ exploits 27 P-DIS data sets [NMC, E139, EMC], 6 p+A DY data sets [E772, E866], together with 6 v-DIS 
data sets [NuTeV, CDHSW, CHORUS] and 3 pion data sets in d+Au at RHIC [PHENIX,STAR], summing up to 
1579 datapoints altogether (for the Refs., see Il20l ). The DSSZ analysis thus covers the most extensive selection of 
the nuclear hard-process data types so far. The key features, and simultaneously also the key differences relative to 
EPS09, are that no weights are introduced for the data sets, and that in computing the inclusive pion production in 
d+Au collisions at RHIC, DSSZ uses the nuclear fragmentation functions (nFFs) determined in [36]. The effects of 
uncertainties in the free-proton MSTW PDFs are not examined in nuclear ratios but DSSZ cleverly includes these 
as additional errors to the measured absolute v-DIS structure functions. The DSSZ analysis arrives at a very good 
overall fit, ^- 2 /datapoint = 0.98. The propagation of the experimental uncertainties into the nPDFs is quantified via 
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the Hessian analysis, which yields 25+25 error sets corresponding to a choice Ay 2 = 30 (i.e. allowing for a 2% 
increase ofx 1 )- The conclusion from the DSSZ analysis is that good simultaneous and tensionless fits to the data used 
are found, and, in particular, there seems to be no conflict between the nuclear modifications of the / ± -DIS and v-DIS 
data. Also based on DSSZ, the measured nuclear effects in hard processes can be explained by universal nPDFs. For a 
detailed account of the DSSZ nuclear modifications and their error bands, see Figs. 13-17 in Ref. |20|. A comparison 
of the DSSZ with the EPS09 and nDS nuclear effects at a scale Q 2 = 10 GeV 2 is presented in Fig. 14 of Ref. EDI . 




10" 4 10 J 10' 2 10"' 1 10 J ID' 2 10 J 1 10 J 10 J 10"' 1 

XXX 

Figure 3: (Color on-line) Upper panels: The DSSZ and EPS09 nuclear effects to bound-proton PDFs (valence u quarks, sea u quarks, gluons) in a 
lead nucleus as a function of x at the initial scales 2q = 1.0 GeV 2 for DSSZ (blue long dashed) and 1 .69 GeV 2 for EPS09 (green solid). The error 
bands are from EPS09. The DSSZ modifications at Q 2 = 1.69 GeV 2 are also shown (red short dashed); their complicated behaviour is caused by 
the negative gluon PDFs in MSTW and DSSZ at smallest x and small Q 2 . Lower panels: The same but for Q 2 = 10 GeV 2 . Prepared by I. Helenius. 

It is interesting to compare the EPS09 and DSSZ results, see Fig.[3]above, especially in the gluon shadowing and 
antishadowing regions, where the results from these two analyses differ the most. The initial nuclear modifications 
Rt(x, Qq) in the smallest-* region which is not directly constrained by the existing nuclear hard-process data, are 
fairly similar in DSSZ (long dashed) and in EPS09 (solid) but their different initial scales should be noticed. The 
negative smallest- x gluon PDFs at low scales both in MSTW and in DSSZ (which as such are not problematic, the 
PDFs themselves evolve in a stable manner and the physical cross-sections remain of course positive) combined with 
the rapid DGLAP evolution of the small-x gluons at Q ~ 1 GeV obscure the interpretation of the R^(x, Q 2 ) in DSSZ at 
the smallest scales above (2o- As shown in the upper right panel on the right in Fig.[T]illustrating these scale evolution 
effects, with the non-negative gluons in EPS09 and CTEQ6M these interpretation (and control) problems do not occur. 

As seen in Fig. [3] in the valence and sea quark sectors the DSSZ and EPS09 do not differ significantly. At the 
smallest x, the sea-quark differences reflect the lack of small-x, high-Q 2 nuclear DIS data. For the valence quarks, 
there is a difference in the mid-x EMC -region, the origin of which should be clarified. 

Figure|3]also shows that the amounts of gluon antishadowing at x ~ 0. 1 in DSSZ and EPS09 are quite different, and 
while the EPS09 gluons show a clear EMC effect there is none in DSSZ. The EMC effect and strong antishadowing 
in the EPS09 gluons follow essentially from the PHENIX n° minimum-bias data OTI (see Fig.|2]above), which show 
these effects and which were given a fairly significant weight in the lack of other gluonic constraints in this region. As 
far as I can see, the explanation for the lack of these effects in the DSSZ gluons is that they are already fully included 
in the nFFs [36| which DSSZ uses to compute the pion production cross sections. What seems troublesome to me 
here is that in the nFF analysis [36 1 and in DSSZ one uses essentially the same pion production data to fix the nuclear 
effects in the gluonic nFFs. Since the nDS nPDFs, which do not show gluon antishadowing or EMC effects either, are 
utilized in the determination of these nFFs, essentially all the antishadowing and the beginning of the EMC effect in the 
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measured pion production are transferred into modifications of the gluonic nFFs. Then, by construction, these effects 
cannot anymore arise for the nPDFs in the DSSZ analysis, either. Clearly, this situation calls for a combined global 
analysis for the nFFs and nPDFs, which would exploit the different types of inclusive DIS, DY and pion production 
data together with the SIDIS data (see the Refs. in (36)). In addition, one should also consider whether the very large 
nuclear effects obtained in the nFFs of Ref. [36] could be understood in terms of some feasible physical mechanism. 

1.5 

1.4 
1.3 
1.2 

%| 1.0 
0.9 
0.8 
0.7 
0.6 

2 4 6 8 10 12 14 16 18 20 
p T [GeV/c] 

Figure 4: (Color on-line) The minimum-bias nuclear modification ratio fij^ u as a function of pj, see the text for details. Prepared by I. Helenius. 

Figure [4] illustrates the rather significant role the nFFs play in the DSSZ analysis. The PHENIX minimum-bias 
data for the nuclear modification factor R^ A for 7r° production in d+Au collisions at RHIC, is scaled up by the factor 
1.039 as in the EPS09 analysis (and also in Fig.|2|l. The EPS09 result (green solid curve and error band) computed 
here with the CTEQ6M PDFs and fDSS vacuum FFs [37|, by construction fits the data. The corresponding result 
obtained with the DSSZ nPDFs and nFFS, using the CTEQ6M PDFs, is shown by the red dashed, and the result 
for DSSZ+nFF with MSTW by the black dotted curve (the p+p baseline in all these cases was computed with fDSS 
FFs). Unlike the nFFS, the free-proton PDFs are thus not a source of a large uncertainty here. The effect of changing 
the vacuum FFs into nFFs is, however, large: with DSSZ+vacuum FFs (blue dashed, purple dotted) essentially all 
the antishadowing and EMC effects in the computed tt° production are now gone, indicating that all these effects are 
essentially included only in the (especially gluonic) nFFs in the DSSZ analysis. 
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5. Impact-parameter dependent nPDFs: EPS09s and EKS98s 

The spatial (impact parameter) dependence arises naturally in the modeling of the origin for the nuclear shadowing, 
see e.g. E51 . but due to the lack of sufficient data constraints it has so far not been possible to embed such a dependence 
directly into the global analysis of nPDFs. However, for the computation of hard-process cross sections in different 
centrality classes (impact parameters b) of nuclear collisions, 

dN AB^ k+ x (b) = £_L\- fd 2 Sl T A ( Sl ),f(xuQ 2 , S i)f^(xuQ 2 ) 

i,j,X' N A ,N B J 

® j dh 2 T B (s 2 )^( X2 M\s 2 )f^(x 2 ,Q 2 )®d& i ^ k+x '6 (2) (s 2 -s 1 -b), (3) 

one quite clearly needs to know also the spatial dependence of the nPDFs, i.e. how different the nPDFs are near 
the edge of the nucleus (s ~ Ra) from those at the center (s ~ 0). In the expression above, we have introduced 
the spatially dependent nuclear modification rf(x, Q 2 ,s) of the PDFs which is related to the average modification 
Rf(x, Q 2 ) for each parton flavor i by 

Rf(x, Q 2 )=jj d2 s T A (s) rf(x, Q 2 , s), (4) 
where Ta(s) is the standard nuclear thickness function. The indices Na(Nb) run through the nucleons in A(B). 



Kari J. Eskola /Nuclear Physics A 00 (2012 ) 74S] 



7 



As discussed by Helenius in these proceedings [22], further light on this issue was recently shed by the study 
in Ref. GTTl . where the A dependence of the EPS09 and EKS98 globally analysed nPDFs was turned into a spatial 
dependence through an ansatz involving powers of Ta(s), 

n 

rf{x, Q\ s) = 1 + J] C>. G 2 ) [Ta(s)V . (5) 

7=1 

The main idea here is that the A-independent coefficients c'j(x, Q 2 ) are obtained for each x, each Q 2 , each flavor i of 
a particular nPDF set by performing a fit to reproduce the normalization in Eq. |4]) as a function of A. As shown in 
ETII the (measured) A systematics of the nPDFs can be reproduced nicely with n = 4. For figures illustrating the 
resulting spatial dependence of the nPDFs, see Ref. [21 1 and Fig. 2 in [22|. The spatial fits were prepared for the 
EKS98 LO nPDFs, and for the EPS09 LO and NLO nPDFs as well as for all the EPS09 error sets. The outcome of this 
useful exercise, the spatially dependent nPDF sets named "EKS98s" and "EPS09s", is now publicly available^ With 
these, one may now for the first time compute nuclear hard-process cross-sections and estimate their nPDF-originating 
uncertainties in different centrality classes (or impact parameters) consistently with a globally analysed nPDFs. 

As shown in Ref. ll2D and discussed by Helenius in these proceedings 11221 (see Fig. 3 there), with the EPS09s 
nPDFs one can, within the experimental and theoretical uncertainties, nicely reproduce the PHENIX data OTI on 
the centrality dependence of the nuclear modification factor /?J A . Predictions for the centrality dependence of the 
corresponding modification in the forthcoming p+Pb collisions at V*aw = 5.0 TeV at the LHC have also been 
computed with EPS09s, see Ref. ETI . and Fig. 4 in El . 

6. Conclusions 

Based on the recent global analyses of nPDFs, I conclude the following: f i) The very good quality fits obtained 
suggest that in the x, Q 2 -region probed by the data, x > 0.001 and Q 2 > 1 GeV 2 , factorization works and the 
nPDFs seem universal, (ii) The NLO analyses with error sets (EPS09, DSSZ) have brought the nPDF global fits to 
the similar (NLO) sophistication level as the free-proton PDF analyses. (Hi) There are still large uncertainties in the 
gluon sector: The role of the nuclear FFs in understanding R^ Au should be clarified by performing a simultaneous 
global fit of nPDFs and nFFs. (iv) The free-proton PDF uncertainties are so far only partially (or not at all) accounted 
for in the nPDF analyses. Ultimately, one should develop a combined global analysis of PDFs and nPDFs. (v) 
Further hard-process data from RHIC d+Au and LHC p+Pb collisions, such as direct photons, high-pj pions, heavy 
quark+photon production and Z/W asymmetries, will help in constraining the nPDFs further. ( vi) To resolve the gluon 
uncertainties at small-x&high-Q 2 , DIS data from future e+A colliders (EIC, LHeC) would be needed. ( vii) The impact 
parameter dependence of globally analysed nPDFs has now been determined (in a specific framework) for the first 
time consistently with their A-systematics: with EPS09s (NLO&LO + error sets) and EKS98s (LO) it is now possible 
to compute nuclear hard-process cross-sections in different centrality classes consistently with EPS09 and EKS98. A 
future task is to implement this type of framework directly into the global fits. 

Acknowledgments. I thank Ilkka Helenius for preparing the Figs. [3] and [4] above, T. Lappi and I. Helenius for 
comments regarding the talk and this manuscript, and DSSZ for providing us with the DSSZ nPDFs. The financial 
support from the Academy of Finland, project 133005, and Magnus Ehrnrooth Foundation is gratefully acknowledged. 
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